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1. Introduction 

The TAUOLA package |1I2I3I4| for the simu- 
lation of r-lepton decays and PHOTOS |5I6| for 
the simulation of radiative corrections in decays, 
are computing projects with a rather long his- 
tory. Written and maintained by well-defined 
authors, they nonetheless migrated into a wide 
range of applications where they became ingredi- 
ents of complicated simulation chains. As a con- 
sequence, a large number of different versions are 
presently in use. From the algorithmic point of 
view, they often differ only in a few small details, 
but incorporate many specific results from dis- 
tinct r-lepton measurements or phenomenological 
projects. Such versions were mainly maintained 
(and will remain so) by the experiments taking 
precision data on r leptons. On the other hand, 
many new applications were developed recently, 
often requiring a program interface to other pack- 
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ages (e.g. generating events for LHC, LC, Belle 
or BaBar physics processes). The programs struc- 
ture, prepared for the convenience of users, was 
presented during previous r conferences, and we 
will not repeat it here. 

This time, let us concentrate on physics ori- 
ented results and special techniques being devel- 
opped and used for fits, over the years. They 
may be worth some attention now, when high 
precision, high statistic data are being analyzed. 
FORTRAN to C++ migration is not as simple as one 
may think. On one hand it is generally believed 
that large projects are easier to maintain in more 
modern programming frameworks. On the other 
hand, FORTRAN software can not be abandoned. 
It usually carries a lot of physics expertise. It 
is also impractical to replace the software envi- 
ronment of mature experiments with something 
new. Fortunately, co-existence of the two lan- 
guages is not a problem, at least not from the 
software point of view. 

Our presentation is organized as follows. Sec- 
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tion 2 is devoted to the discussion of optional 
weights in TAUOLA and their use for fits to the 
data already at the level of comparison with 
raw data. In section 3 we present some new 
results for the simulation with PHOTOS of ra- 
diative corrections in decays. Section 4 is de- 
voted to MC-TESTER; the program which can be 
used for semi-automatic comparisons of simula- 
tion samples originating from different programs. 
The prototype, C++ based, TAUOLA universal 
interface is mentioned in that Section too. 

Because of the limited space of the contribu- 
tion, and sizable amount of other physically inter- 
esting results, some of them will be excluded from 
conference proceedings. Let us hope that they 
will find place in future works, possibly with col- 
laborators mentioned in the Abstract. For these 
works, the present paper may serve as an an- 
nouncement. 

2. Optional weights in TAUOLA Monte Carlo 

Physics of r lepton decays requires sophis- 
ticated strategies for the confrontation of phe- 
nomenological models with experimental data. 
On one side high statistics experimental samples 
are collected, and the obtained precision is high, 
on the other hand, there is a significant cross- 
contamination between distinct r decay channels. 
Starting from a certain precision level all chan- 
nels need to be analyzed simultaneously. Change 
of parametrization for one channel contributiong 
to the background may be important for the fit 
of another one. This situation leads to a com- 
plex configuatin where a multitude of parame- 
ters needs to be simultaneously confronted with 
a multitude of observables. One has to keep in 
mind that the models used to obtain distribu- 
tionbs in the fits may require refinements or even 
substantial rebuilds. The isospin relation to dis- 
tribution measured in e + e~ collistions at low en- 
ergies offer additional control, but may be the 
source of even more complex systematic error con- 
siderations. Effects of QED bremsstrahlung need 
to be taken ito account as well. It may affect 
selection criteria and background contaminations 
in quite complex and unexpected ways. 

From the statistical point of view it is best to 



resolve such system in one automated step using 
the method such as |7l8j for example. This can be 
of course very dangerous from the point of view of 
systematic error control. We will return to this 
point later. At this moment let us discuss nec- 
essary adaptations for TAUOLA for such methods 
to be used. In fact the principle of the changes 
is quite simple. The method is in use since long 
time and in different situations, see eg. [9]. It is 
enough, to calculate, for each generated event (for 
each present in it decay of r + and t~ separately) 
alternative weights; the ratios of matrix element 
squared obtained with new currents provided by 
the user, and the one actually used in generation 
are then the vector of weights. If the currents can 
be represented as a linear combination of elemen- 
tary ones, then the weight for arbitary values of 
coefficients is a bilinear form in the relative coef- 
ficient^- In more general cases, linearization is 
necessary, and the whole procedure needs to be 
repeated with the currents obtained at one step 
used for generation at the next one. In this way, 
effects of the new currents will migrate to spin 
correlations between r + and t~ as well. 

The first example of the use of the method, is 
demonstrated in fig. [1] At this step of the work 
no detector effects were included. Experimen- 
tal data were provided as a one-dimensional his- 
togram. Simulation of r production process was 
not performed. One can see, that the theoretical 
model provides double peak which is absent in the 
data. This is because of an interplay between in- 
termediate state resonance parametrization and 
phase space jacobian. The same parametrization 
was used for 3-scalar and 2-scalar systems. This 
points to the necessity of further investigations 
of unitarity constraints on the running width or 
better control of the data. 

Options how to handle, in experimental soft- 
ware, information for weight calculation is under 
study now. In principle, information stored in 
the event record should be enough. This may 
not be a perfect match in all cases. For exam- 
ple, effects of the boost from laboratory frame to 
rest-frame of r lepton and back, may result in 

1 ~We have applied such method in one of our other pro- 
grams [9] for the studies of anomalous couplings of Z and 
LEP data. 
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Figure 1. Preliminary result from the study of t~ — > 
rj-K^-K ^ decay using external weight. Experimental 
data and theoretical predictions obtained at the inter- 
mediate step of the work are shown. The distribu- 
tion of the invariant mass of all final state hadrons 
is given. This plot is available thanks to courtesy of 
Vladimir Cherepanov. 



rounding errors. In some cases the kinematical 
configuration may be changed because of QED 
bremsstrahlung generated by PHOTOS. The order 
of final state particles stored in the event record 
may not be the same as used internally in TAUOLA 
at the event construction step. This of course can 
be cross checked on the basis of inspection of a 
smaller sample and without detector simulation 
part being used. Such tests are nonetheless nec- 
essary, but may need less time than would be the 
case if collaboration software would be rebuild. 
From that point of view, the easiest would be to 
calculate optional weights from information al- 
ready available in experiments production tapes, 
and only for those events where bremsstrahlung 
effects are not generated at all. This can be cor- 
reted at the next step of form-factor iteration. 
Related slow down of the convergence should not 
be large. Possibly two steps will be enough, the 
second one for the verification. It will be neces- 
sary for control of spin correlations between the 
two decaying r-s as well. It is technically much 



simpler not to include changes into spin corre- 
lations. In such a case the vector of weights 
would depend on kinematical configurations of 
r + , r~ and of the production process simulta- 
neously. The optimal choice require further in- 
teraction with collaborations. Technical side of 
the activity will depend on the way the user want 
to calculate the hadronic current. It can be easily 
calculated as external generalized function com- 
municating with TAUOLA at the level of linker or 
through the f if o pipes. They can be provided as 
FORTRAN or eg. C++ code and have well defined 
and easy to control input/output for modifica- 
tions introduced into TAUOLA. 

Int the past, and for a good reason, signifi- 
cant effort was devoted to strategies where the 
methotds as indicated here can be avoided. For 
example in [lOj a system of projections was pro- 
posed for the 3-scalar r decays. If such a method 
can be used it is possible to obtain from the data 
just one of the contribution to the current, which 
is then parametrized by a single scalar function of 
two invariants only. This is of course on the ex- 
pense of statistical strength on the measurement 
and abstraction of experimental sytematic error. 
Such methods, easy to apply for final states of 
two or three scalars, are less practical for higher 
multiplicities. Nonetheless they can be of great 
help to control the analysis and to verify which 
aspects of the models require improvements. 

3. PHOTOS and NLO effects in W and 7* 
7r + 7r decays 

The changes introduced over the last two years 
into PHOTOS Monte Carlo program itself were 
rather small. On the other hand, the complexity 
of the work on its theoretical foundation was large 
and matches neither size nor the main purpose of 
the present talk (some of the related topics are re- 
vieved in [H]). Previous tests of two body decays 
of the Z into a pair of charged leptons [H] and 
scalar B into a pair of scalars [13] was recently 
supplemented [14] with the study of W r± — > l^vf. 
The study of the process for 7* — > it + it~ is on- 
going 15J. In all of these cases universal kernel of 
PHOTOS was replaced with the one matching exact 
first order matrix element. In this way terms for 
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NLO/NLL level are implemented. As usual, algo- 
rithm covers the full multiphoton phase-space and 
is exact in the infrared region of the phase space. 
This is rather unusual for the NLL compatible al- 
gorithms. One should not forget that PHOTOS gen- 
erates weight one events and does not require any 
factorization scales. There is full overlap of the 
phase space where hard matrix element is used 
with the one for iterative photon emissions. All 
interference effects are implemented with the help 
of internal weights. 

The results of all tests of PHOTOS with a NLO 
kernel are at sub-permille level. No differences 
with benchmarks can be observed, even for sam- 
ples of 10 9 events. For the comparisons when 
simpler physics assumptions are used, differences 
were sub-permille too, if compared to total rates 
(or were matching precision of programs used for 
tests). 

This is very encouraging, and points to the 
possible extension of the approach outside of the 
QED (scalar QED) only. In particular to the do- 
main of QCD or if phenomenological lagrangians 
for interactions of photons need to be applied. 
For that work to be completed spin amplitudes 
need to be further studied. Let us point here to 
the paper [16] . 

The discussed here refinements only indirectly 
affect practical side of simulations for r physics. 
Changes in the kernels necessary for NLO may 
remain as options for tests only. They are avail- 
able from the PHOTOS web page [14] but are not 
recomended for wider use. An example of such 
new tests is given in fig. Eb- 

4. MC-TESTER and C++ TAUOLA universal 
interface 

In the development of packages such as TAUOLA 
or PHOTOS, question of tests and appriopriate re- 
lations to users applications are essential for their 
usefulness. In fact, users applications may be 
much larger in size and human effort than the 
programs discussed here. Good example of such 
'user applications' are complete environments to 
fit experimental data and control detector uncer- 
tainties at the same time. As a consequence our 
programs need to work well with FORTRAN where 



HEPEVT event record is used and C++ where HepMC 
[17] is used instead. For the sake of automatiza- 
tion of tests MC-TESTER was prepared [18]. Re- 
cently, new functionalities were introduced into 
the package 19J. In particular, it works now 
with the HepMC event record, the standard of 
C++ programs. The complete set-up for bench- 
marking the interfaces, such as interface between 
r-lepton production and decay, including QED 
bremsstrahlung effects is prepared. The exam- 
ple is chosen to illustrate the new options intro- 
duced into the program and novel ways of its use. 
From the technical perspective, [19] documents 
program updates and supplements previous doc- 
umentation [T5] , 

As in the past, test consists of two steps. Dis- 
tinct Monte Carlo programs are run separately; 
events are searched for the decays of the chosen 
object, and information is stored by MC-TESTER. 
Then, at the analysis step, information from a 
pair of such runs may be compared and repre- 
sented in the form of tables and plots. 

One may have the impression that the modifi- 
cations introduced into the release are minor and 
consist of simple improvements in the graphical 
representation of the output and purely technical 
reorganization thanks to the use of C++ lists. 

To some degree this is true, but on the other 
hand the changes were introduced because of 
pressure from applications. Let us show how pro- 
gram modifications can be used for non-trivial 
practical applications. 

It is quite common that information stored in 
the event record is too large. For example in- 
dividual soft photons which remain undetectable 
are present. Not only they do not influence the 
detector response at all, but they exhibit tech- 
nical aspects of eg. infrared regulators of QED 
bremsstrahlung. In response, MC-TESTER should 
ignore those photons while analyzing decays (or 
group together with other particles). Other- 
wise, comparisons of different Monte Carlo pro- 
grams would be dominated by the technical as- 
pect of the implementation of the infrared regu- 
lator; MC-TESTER operation need to be adopted 
for this, see eg. |20|12| . 

Another example where the event tree may 
need to be simplified for validation is if spin corre- 
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Figure 2. Examples of interesting benchmarks: 
In the decay Z — > t + t - — > tx^-k^vv, the tt~ en- 
ergy spectrum in the Z rest-frame has identical 
shape as the distribution of the invariant mass 
squared ofTr + ir~D, see plot a. This linear distri- 
bution can be used to measure decaying t~ polar- 
ization. In plot b, an example of the PHOTOS test 
is given. Program version with full matrix ele- 
ment for W + — > /i + i^,7 channel is compared with 
the standard one. The distribution of the /i + 7 
pair mass is shown. 



lations are appropriately introduced into various 



production processes. Let us use as an exampld3 
pp —> Z/"f * +X, Z — » t + t~. It is convenient to 
start the test by restricting t decays to the sim- 



plest decay mode, that is r d 



and look 



at distributions in cascade decay Z — > ■k + tt~vv. 
In this case the effects of spin correlations are 
largest. The distribution of the 7r _ energy spec- 
trum (in the Z rest- frame), manifests the r po- 
larization through its slope (see fig [5^) . Fortu- 
nately, this frame dependent spectrum is equal to 
the distribution of the square of invariant mass 
of Tt~ir + i>. This distribution can be obtained in 
MC-TESTER thanks to the new options. 

Final state activities will lead to final states 
of Z decays where the r-pair is accompanied by 
bremsstrahlung photons or soft hadrons. Such 
distribution are usually strongly peaked. In 
fig. Eb we provide an example of such a distribu- 
tion but for the comparison of PHOTOS with exact 
marix element for the decay W + —> [i + v^f. 

With the above brief examples we have demon- 
strated that MC-TESTER may be useful for tests of 
libraries of particles decays, as well as for tests 
of their interfaces. The work on HepMC based 
TAUOLA universal interface is on-going. A 
preliminary version is already available [21j . but 
from physics point of view it is still less complete 
than the FORTRAN version [22]. 

The updated version of MC-TESTER was found 
|20|12| to handle well cases where physically spu- 
rious information (eg. on soft photons) need to 
be ignored. 

5. Summary and future possibilities 

The status of the computer programs for the 
decay of r leptons TAUOLA and associated projects 
TAUOLA universal interface and MC-TESTER 
was reviewed. The high-precision version of 
PHOTOS for radiative corrections was presented 
too. 

New results for PHOTOS and decay of W were 
mentioned. For this channel complete next-to- 
leading order effects can now be simulated. How- 
ever, for most of the applications these effects 
are not necessary, leaving the standard modu- 



2 Spin correlations in decays of W, H, into r lepton(s) 
are nearly identical. 
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lar version of PHOTOS sufficient. This exercise is 
important not only for photonic bremsstrahlung. 
It helps to understand better questions related 
with matching hard emission matrix elements 
with parton showers without necessity to intro- 
duce any boundaries within phase space. The 
important result of the above work is that the 
path to include electromagnetic form-factors of 
the particles participating in decay is now open 
for future fits to the data. These form-factor ef- 
fects may be significantly larger and physically 
more justified than complete next-to- leading or- 
der effects of scalar QED as in B meson decays. 

The presentation of the TAUOLA general- 
purpose interface was brief. The C++ version ex- 
ist, but is still physically less refined than the 
FORTRAN version. Examples for its use can be 
found on the continuously modified web page [21] . 
The new version of MC-TESTER [19] is already pub- 
lic now. It not only works with HepMC [T7] of C++ 
but enables user defined tests. 

Distinct versions of the TAUOLA library for t 
lepton decay, and of PHOTOS for radiative correc- 
tions in decays, are now in use. The principles 
how to use the distribution package are presented 
in refs. [4123] , on the other hand, development of 
new and recovery of codes for old strategies of 
fits is becoming available. This includes strategy 
to use weighted events, and the projection oper- 
ators [10] for the three scalars final states. Let 
us hope that the present talk will provide some 
stimulation for the larger work in the near future. 
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